Electron transport through a diatomic molecule 
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Abstract 

Electron transport through a diatomic molecular tunnel junction shows wave like interference phenomenon. By using Keldysh 
non-equilibrium Green's function (NEGF) theory, we have explicitly presented current and differential conductance calculation 
for a diatomic molecular and two isolated atoms (two atoms having zero hybridization between their energy orbital) tunnel 
junctions. In case of a diatomic molecular tunnel junction, Green's function propagators entering into current and differ- 
ential conductance formula interfere constructively for a molecular anti-bonding state and destructively for bonding state. 
Consequently, conductance through a molecular bonding state is suppressed, and to conserve current, conductance through 
anti-bonding state is enhanced. Therefore, current steps and differential conductance peaks amplitude show asymmetric corre- 
spondence between molecular bonding and anti-bonding states. Interestingly, for a diatomic molecule, comprising of two atoms 
of same energy level, these propagators interfere completely destructively for molecular bonding state and constructively for 
molecular anti-bonding state. Hence under such condition, a single step or a single peak is shown up in current versus voltage 
or differential conductance versus voltage studies. 

PACS numbers: 73.23. Hk, 85.85. 85.65.-|-h 
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I. INTRODUCTION 

Aptitude of semiconductor technology towards minia- 
turisation of electronic devices has lead to emergence of 
nano and molecular electronics. Now the work of nano 
scientists is to control electron transport at such minute 
scale. As the size of the nano devices is comparable to the 
wavelength of the current carriers so quantum features 
are playing dominant role in these devices. A theoreti- 
cal model of electron transport through a single molecule 
was first presented by A.Aviram and Mark.A.Ratneri'^. 
Since then there has been great interest both from exper- 
imental and theoretical groups for the progress of molec- 
ular electronics^"—. A single molecule could be used as 
electronic mixersS, switches^*^, and rectifiers^. With an 
eye on its applications, it is expected that the under- 
standing of quantum electron transport at the molecular 
scale is a key step towards practical usage of molecular 
devicesi^. Experiments on conduction through molecular 
tunnel junctions are becoming more common, for discus- 
sion see Reiiiii^ and references therein. On the experi- 
mental front, the most common methods of contacting in- 
dividual molecule are scanning tunnelling microscope tip 
and mechanically controlled break junctionsi^ii^. Early 
experiments in the field of molecular electronics have fo- 
cused on the absolute conductance and its dependence 
on wire length, molecular structure, and temperature-^. 
From the theoretical point of view, investigating electron 
transport in an electrically contacted molecule is a chal- 
lenging problem. Most of the formal theoretical work on 
transport through molecular tunnel junction have relied 
on the Generalized master equation's theorjsi^ii^ and the 
non equilibrium Green's function (NEGF) theoryi^ii^. 
By using NEGF theory a plethora of physical phenomena 
like coulomb blockadei^, Kondo effect^ and vibrations 
effeclj^ are already reported. 

Electron transport through a molecular tunnel junc- 
tion has very non-intuitive characteristics. Many experi- 
mentalists and theoreticians of molecular science refer to 



hydrogen molecular transport to be explained by trans- 
port through a single channel^"—, instead of the pres- 
ence of hydrogen molecular bonding and anti-bonding 
states, as conductance channels. On the basis of the first- 
principle calculations, it has been already reported that 
the transmission probability through hydrogen molecular 
tunnel junction becomes exactly one for molecular anti- 
bonding stato^. Apart from the hydrogen molecule, even 
more complex molecules show asymmetry in conductance 
through highest occupied molecular orbital (HOMO) and 
lowest unoccupied molecular orbital (LUMO)^'^. 

In this study we discuss a diatomic molecular tunnel 
junction and contrast it with two isolated atoms tunnel 
junction. The results obtained by the present theoretical 
formulation is valid even for artificial diatomic molecule 
(two coupled quantum dots'^'^). By using NEGF theory, 
we have explicitly shown the current and differential con- 
ductance of the under discussion system. In the presence 
of the hybridization between energy orbital of the two 
atoms, a diatomic molecular bonding and anti-bonding 
states are formed. The Green's function propagators 
entering into current and differential conductance inter- 
fere constructively for two isolated atoms tunnel junc- 
tion, where as, for a diatomic molecular tunnel junc- 
tion, they interfere constructively for an anti-bonding 
state and destructively for a bonding state. There- 
fore, the normalized transmission probability function 
for a diatomic molecular tunnel junction shows asym- 
metric correspondence between molecular bonding and 
anti-bonding states. For a diatomic molecule compris- 
ing of two atoms of the same energy level the normal- 
ized transmission probability becomes one for molecular 
anti-bonding state and zero for molecular bonding state. 
This clarifies the asymmetric correspondence in current 
steps and differential conductance peaks amplitude for a 
diatomic molecular tunnel junction. 

This paper is organized as follows. In section II a 
diatomic molecular tunnel junction model Hamiltonian, 
current and differential conductance are presented. In 
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section III numerical results are discussed. Conclusions 
are made in section IV. 



II. THEORETICAL MODELLING 



of the thermal energy, with Boltzmann's constant ks and 
temperature T. By employing wide band approximation 
{J2k — ^ ^ I^oo ^^'^ using electrode Green's func- 
tion in Eq. (3), the current from a diatomic molecular 
tunnel junction is given by the following relation. 



In modelling Hamiltonian of a diatomic molecule em- 
ployed between electrodes, we have considered electrodes 
as charge carriers bath and charge carriers wave function 
hybridize between electrodes and i-th atom of a diatomic 
molecule. The Hamiltonian of the under discussion meso- 
scopic system is, 

2 2 

H = Y^ eidjd, + Td\dj + 

kui k 

The Hamiltonian is written under the approximation of 
linear combination of atomic orbitals and in second quan- 
tized form^^. Here dl{di) represents creation (annihila- 
tion) operator of electron on i-th atom, where as second 
term refers to inter-atomic hybridization energy r. The 
CijC^k, and T^k are representing i-th atom, electrodes, and 
i-th atom and electrodes hybridization energy, respec- 
tively. The operator cli.{ci^k) is electron creation (anni- 
hilation) operator on the ly electrode. We have assumed 
for simplicity h = e = 1. 

The current across a diatomic molecular tunnel junc- 
tion is calculated by the time evolution of the occupation 
number operator n"^ (t) = (t) cli^c^k of the v elec- 

trode. 



= 2Re^T,krGc<{t,t)] 
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Here, Gc^{t,t) — i {^c^j^dij represents electrode and a 

diatomic molecular coupled lesser Green's function. A di- 
atomic molecular and v electrode coupled lesser Green's 
function can be de-coupled into a diatomic molecular 
(G^'^) and v electrode (ff^^^) Green's function by ex- 
ploiting analytic continuation rul( 



,18 



juk 



(3) 



+G<{t,h)g^,{h,t')] 



The v electrode lesser and advanced Green's function 
are g^^{t, t') = i6 {f - t) Exp [ze^ {t - t')] and g<.(t, t') = 
ifu (cfc) Exp [~iek {t - t')]. 

Here (e^) = — r— rr is the Fermi-Dirac 

Exp [/3 (ck - Ml + 1 

distribution function. Where as, /? — 
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gives inverse 



/„ (t) = -2 In 



(Exp [ick (t - ti)]) [/„ (efc) Gg(t, ti) -I- G< (t, ti)]] 



(4) 



Here T^ij — 2-KDTukiTl^y By utilizing the conditions of 
current continuity and line width proportionality^ ^Lij — 
V^Riji the current through a diatomic molecular tunnel 
junction is given as follows. 



Il (t) = 



Ul (efe) - Sr Exp [iefc (t - ti)] Gfj(t, ti)] 



(5) 



Here, T = Tmj = Tnj. Now we use Dyson's equation to 
find a diatomic molecular-system total retarded Green's 
function Gf^. 



G5(t,t') =ff5(t,t') + J^T„,.Tt, 



I' dti J dt2s5(t,ti)gffc(ti,t2)Gfj(t2,t') 



(6) 



Here, and g^f, represent a diatomic molecular and f 
electrode retarded Green's functions, respectively. By 
using electrode Green's function into Eq.(6), diatomic 
molecular-system total retarded Green's function is sim- 
plified. 



GM,t')^gm,t')Exp[-T{t-t')] 



(7) 



(2) Now a diatomic molecular retarded Green's g^ function 
is found by the equation of motion technique. 

d 

i—[Exp{ie,t)di{t)] =TExp{ie,t)d2{t) (8) 



i — [Exp{i€2t)d2{t)] = TExp{ie2t)di{t) (9) 

By defining Fi{t) = Exp{ie^t)di{t) 
and F2{t) = Exp{ie2t)d2{t) 

-Exp{iM)-E{t) ^ F2{t) (10) 
r ot 



-Exp{~iM)—F2{t) = Flit) (11) 

T ot 

Here A — {e^ — e^). For de-coupling Eqs. 10 and 11, we 

operate ?— on above two equations. 
ot 



(12) 
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2 . ^ dF2 
dt^ dt 



T^F2 = 



(13) 



By solving above two differential equations, we find time 
evolution of the annihilation operators of the diatomic 
molecule. 



diit) = — [{(7 + A)di(O) -2rd2(0)}£;xp(-ie6t)+} 
27 

{(7 - A)di(O) + 2Td2(0)}Exp{-ieat)}] 



d2{t) = —[{(7 - A)d2{0) - 2Tdi{0)}Exp{~iett)+ 
27 

{(7 + A)(i2(0) + 2Tdi{0)}Exp{~ieat)}] 

(ei + £2 - 7) 



(14) 

(15) 
and 



Here 7 = y/A^ + ir^ , e& 

z 

Ea = + ^2 + 7) ^ Where as e6(ea) represents molec- 
ular bonding (molecular anti-bonding) energy. Similarly, 

= T7-[«'^ + ^)4(0) - 2Tdl{0)}Exp{ie,t)+ 
27 

{(7 - A)4(0) + 2Tdl{0)}Exp{ieat)}] 



(16) 
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^[{(7 - A)4(0) - 2Tdl{Q)}Exp{ie,t)+ 
27 



{(7 A)4(0) 2r4 (0)}£;xp(i£ai)}] (17) 

Now a diatomic molecular retarded Green's function gfj 
is found by the following relation. 

g^^it,t')^-^e{t-t'){m) d]it')]) (18) 

By utilizing Eqs. 14-18 and using anti-commutation 
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a diatomic molecular-system 



total retarded Green's function is found. 



Gfa (6) 
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-ea + ir) 


(e - £6 + ir) 


1 r 


(7 - A) 


(7 + A) 



27 ^ (e - £6 + ir) {e-ea+ iF) ' 
Finally, current and differential conductance for a di- 
atomic molecular tunnel junction are given by the fol- 
lowing relations. 



1 = I ^T{e,)[fLie,)^fRie,)] 



(19) 



and 



dV J 2n ^ 



£;xp[/3(€fc - ml)] 
{Exp[(3 {e^ - f,L)] + if 
Exp[l3 (efe - ^ir)] 



{Exp [/3 (efe - pr)] + ly 



with transmission function T (ek) 
r(efc) = -rE.,Im [G«(e,)] 



III. NUMERICAL RESULTS 

The comparative analysis of the electron transport 
between two isolated atoms tunnel junction and a di- 
atomic molecular tunnel junction is discussed in this 
study. The quantum interference effects are usually de- 
stroyed in macroscopic devices. But in nano devices, 
where the wave length of the current carriers become 
comparable with the size of the device, the wave fea- 
tures are preserved. For two isolated atoms tunnel junc- 
tion, current starts flowing as the Fermi energy level of 
the electrode equates with the energy level of the either 
atom. However, when the energy orbital of these atoms 
hybridises then a diatomic molecular bonding and anti- 
bonding states are formed. Now current starts flowing 
as the Fermi energy level of the electrode equates with 
a diatomic molecular bonding or anti-bonding states. In 
this study we have explicitly discussed that the electron 
transport through a diatomic molecular tunnel junction 
shows asymmetric correspondence in current steps and 
differential conductance peaks amplitude for molecular 
bonding and anti-bonding states. More interestingly, for 
a diatomic molecular tunnel junction, where constituents 
of a diatomic molecular energy states are having same 
energy level then only a single step in current and a sin- 
gle peak in differential conductance is shown up^^. This 
could be seen in Fig. 1, where current and differential 
conductance functional relation with applied voltage is 
plotted. When two isolated atoms having the same en- 
ergy level are placed between electrodes then only a sin- 
gle step in current and a single peak in differential con- 
ductance is shown up. While in the presence of finite 
hybridization between energy orbital of these atoms, a 
single step in current and a single peak in differential 
conductance is shown up at the anti-bonding energy level 
of the molecule. Where as, no current step and differ- 
ential conductance peak is shown up at the bonding en- 
ergy level of the molecule. The increase in hybridization 
between energy orbital of these atoms, induces shift be- 
tween molecular bonding and anti-bonding energy levels. 
Consequently, current step and differential conductance 



7 = 1.5 r 1 

r = ? I' ( 

























FIG. 1. (Color online) Current (solid lines) differential con- 
ductance (dashed lines), ei = £2 = 5r and PV — 0.35. 
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FIG. 2. (Color online) Current (solid lines) differential con- 
ductance (dashed lines), ei — SF, £2 = 9T, and /3T = 0.35. 



FIG. 4. (Color online) The normalized transmission proba- 
bility, ei = 5F, and £2 = OF. 



peak too shift to the new molecular anti-bonding energy 
level. 

The asymmetry in conductance through a diatomic 
molecular energy states could be more clearly visu- 
alised by studying electron transport through a diatomic 
molecule comprising of two atoms of different energy 
level. See Fig. 2. In the absence of hybridization be- 
tween energy orbital of the atoms, two symmetric steps 
and peaks are shown up in current and differential con- 
ductance, respectively. The presence of hybridization be- 
tween energy orbital of these atoms results in formation 
of molecular bonding and anti-bonding states. The cur- 
rent step and differential conductance peak amplitude of 
molecular bonding state are suppressed and consequently 
these are enhanced for molecular anti-bonding state. 

To investigate the asymmetry in conductance through 
molecular bonding and anti-bonding states, we discuss 
the normalized transmission probability. The normal- 

ized transmission probability — for two iso- 

lated atoms having the same energy level becomes exactly 
one, when electrodes electron kinetic energy equate with 
energy level of the atom = eo. See Fig. 3. Now 
in the presence of hybridization between energy orbital 
of these atoms, the normalized transmission probability 




-t: 7T maxima shifts to molecular anti-bonding 

energy level Cfe = Eq = to + In the absence of hy- 
bridization between energy orbital of the atoms. Green's 
function propagators interfere constructively and there- 
fore the normalized transmission probability becomes ex- 
actly one, when electrodes electron kinetic energy equate 
with energy level of the atom Ck ~ eo- The Green's 
function propagators for diatomic molecule interfere com- 
pletely constructively for molecular anti-bonding state 
and completely destructively for molecular bonding state. 
This gives perfect transmission through molecular anti- 
bonding state and no transmission through molecular 
bonding state. Therefore, a single step in current and 
a single peak in differential conductance is shown up in 
Fig. 1. 

In Fig. 4 the normalized transmission probability 
for two atoms having different energy level is shown. 
As the kinetic energy of the electrodes electron equate 
with atom of the lower energy level half of the trans- 

. . 1 r2 

vnissinn -_ througli thc devicc is achieved. 



2(efc-ei)2+r2 
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FIG. 3. (Color online) The normalized transmission proba- 
bility. £l = £2 = 5F. 



While the next half of the transmission - , 

2(efc-e2)2+r2 

through the device is achieved where electrodes electron 
kinetic energy equate with atom of the higher energy 
level. In the presence of finite hybridization between en- 
ergy orbitals of these atoms, molecular bonding and anti- 
bonding states are formed. Now the destructive interfer- 
ence of Green's function propagators for molecular bond- 
ing state suppresses the normalized transmission prob- 

. . 1 r2 2r 

ability — — ^(1 — — ), and the constructive 

2 (efc - ea)2 -f 7 

interference of Green's function propagators for molec- 
ular anti-bonding state enhances the normalized trans- 

1 r2 2r 
mission probability ^ t(H )■ Therefore 

2[ek-ear + T^ 7 
current step and differential conductance peak amplitude 
for molecular bonding state is smaller than for molecular 
anti-bonding state, as shown in Fig. 2. With increas- 
ing hybridization between energy orbital of these atoms. 
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the asymmetry in normalized transmission probability, 
current, and differential conductance is increased. 

IV. CONCLUSION 

In this study we have explicitly demonstrated that 
current and differential conductance are suppressed for 
a diatomic molecular bonding state and consequently 
these are enhanced for anti-bonding state. The Green's 
function propagators entering into current and differ- 
ential conductance calculations interfere destructively 
for molecular bonding state and constructively for anti- 



bonding state. Moreover for diatomic molecule compris- 
ing of two atoms of same energy level, no step in current 
and no peak in differential conductance has shown up for 
molecular bonding state. Therefore, for such a diatomic 
molecule current passes through molecular anti-bonding 
state and it does not pass through molecular bonding 
state. And for a diatomic molecule comprising of two 
atoms of different energy level, the current and differen- 
tial conductance have shown to depend upon hybridiza- 
tion between energy orbital of two atoms. The more the 
orbital of the two atoms are hybridised the more current 
passes through molecular anti-bonding state and the less 
current passes through bonding state. 
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